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ABSTRACT 

Observation that DNA molecules in bacteriophage 
capsids preferentially form torus type of knots 
provided a sensitive gauge to evaluate various 
models of DNA arrangement in phage heads. Only 
models resulting in a preponderance of torus knots 
could be considered as close to reality. Recent 
studies revealed that experimentally observed en- 
richment of torus knots can be qualitatively 
reproduced in numerical simulations that include a 
potential inducing nematic arrangement of tightly 
packed DNA molecules within phage capsids. 
Here, we investigate what aspects of the nematic 
arrangement are crucial for inducing formation of 
torus knots. Our results indicate that the effective 
stiffening of DNA by the nematic arrangement not 
only promotes knotting in general but is also the 
decisive factor in promoting formation of DNA 
torus knots in phage capsids. 

INTRODUCTION 

Although structures of proteins forming bacteriophage 
capsids are Icnown with nearly atomic resolution (1-3) the 
precise arrangement of DNA within the capsids remains 
to a large extent unknown. Fluidity of tightly packed 
DNA makes it unsuitable for high-resolution structure de- 
termination based on X-ray crystallography or electron mi- 
croscopy. Although in recent years cryo-electron 
microscopy studies (4,5) and numerical simulations (6,7) 
supported models of coaxial and concentric spooling, 
there is an uncertainty of whether there is a chiral bias in 
formed spools (8). Such a bias could be caused by axial 
rotation of DNA during its active loading into the phage 
capsids (9) or could result from a cholesteric twist angle 



induced by steric interactions between highly concentrated 
right-handed double-helices (10). 

When mature bacteriophages such as P2 or P4 are 
assembled in infected cells their long hnear DNA mol- 
ecules are loaded into phage capsids formed earlier. 
Each phage capsid receives one linear DNA molecule 
that is progressively fed starting with its one end 
through a thin channel that leads through phage tail to 
the phage capsid, which is the main compartment to store 
the phage DNA. The loading process, necessitating 
specialized DNA packaging motors (9), normally stops 
in such a way that the traihng end of loaded DNA is 
retained within the phage tail (4,11), that will serve later 
as the injection needle during the next round of the infec- 
tion, whereas the rest of the long linear DNA molecule 
with its leading end is tightly packed within the phage 
capsid. The two ends of the linear DNA molecules have 
protruding single stranded extensions that are comple- 
mentary to each other (12). This complementarity of 
ends enables rapid circularization of phage DNA upon 
its injection into a newly infected cell. Such a rapid circu- 
larization is crucial for protecting the invading DNA from 
digestion by bacterial exonucleases. However, in tailless 
mutants of bacteriophage P4 the lagging end also enters 
into the phage head during the loading process (13). In 
this specific situation, the complementary ends can 
approach each other as a result of thermal motion 
occurring within the tightly packed but fluid DNA ar- 
rangement in the phage capsid. When the two ends meet 
they anneal with each other so that when the DNA is 
released from the capsid, using standard laboratory pro- 
cedures, the DNA ends remain firmly attached (13). 
Therefore, various types of DNA knots that resulted 
from DNA arrangement in capsids of bacteriophages 
maintain their original knot type after DNA isolation 
and can be investigated by such techniques as electron mi- 
croscopy (13) and gel electrophoresis (14-16). As DNA 
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molecules undergoing cyclization within phage capsids are 
not distorted by any external force, the formed knot types 
can provide information about the overall organization of 
DNA in tightly packed bacteriophage genome. 

Studies investigating relative frequencies with which 
various types of DNA knots formed in bacteriophage 
capsids revealed that these knots have a much different 
spectrum as compared to knots formed randomly. For 
example, while freely fluctuating polymers in free space 
(17,18) or random trajectories highly confined to a 
sphere (8) result about twice more frequently in formation 
of five crossing twist knots (with the standard mathemat- 
ical notation 52) than of five crossing torus knots (5i) (see 
Figure 5 for schematic presentations of these knots) the 
opposite is the case for DNA knots formed in phage 
capsids (8). The standard mathematical notations of 
knots such as 52 use two numbers where the first one, 
written with normal fonts, indicates the minimal number 
of crossings a given knot type can have in a projection 
whereas the subscript number indicates the tabular 
position of a given knot type among the knots with the 
same minimal crossing number in standard tables of 
knots. Another 'anomaly' of knots formed in phage 
heads concerns the four crossing knot 4i (see Figure 5). 
That type of knots due to its relative simplicity forms 
much more frequently than five crossing torus knots (5i) 
or five crossing twist knots (52) by random trajectories 
in confined volumes (8) or by DNA knots formed in 
free solution (19,20). However, among knots formed in 
phage capsids knot 4i is found strongly underrepresented 
and forms significantly less frequently than 5i or 52 
knots (8). 

Although the determination of the specific knot 
spectrum formed by DNA molecules that circularized in 
tailless phages provided a sensitive signature of DNA ar- 
rangement within phage capsids, it is not possible yet to 
directly interpret this knotting signature in terms of pre- 
dicting the geometry of DNA arrangements that will result 
in formation of the same spectrum of DNA knots. The 
researchers in the field rather use a trial and error method 
where they test different models by simulations and 
analyse what type of knots would results from a given 
model. So, for example, Marenduzzo et al. (21,22) have 
shown recently that numerical simulations incorporating a 
potential inducing a nematic arrangement of tightly 
packed DNA molecules produced a very similar spectrum 
of knots to the one observed in real DNA molecules 
that were circularized within bacteriophage capsids. 
Marenduzzo et al. assumed that the cholesteric twist 
angle between sequential layers of nematically arranged 
segments of the same long DNA molecules can be respon- 
sible for the observed preponderance of chiral knots. To 
test this hypothesis, they have compared models where the 
cholesteric twist angle was set to 1°, which is in the upper 
range of the cholesteric twist angle values observed in 
DNA hquid crystals (23), with models where the nematic 
potential was not inducing any cholesteric twist. 
Interestingly, the two conditions produced nearly identical 
spectra of knots. This result indicates that the cholesteric 
twist, which most likely is present in the nematic arrange- 
ment of tightly packed DNA in bacteriophage capsids, is 



not essential for the observed preponderance of torus 
knots in DNA molecules circularized within phage 
capsids. 

What is then the decisive factor in nematic DNA ar- 
rangements that can favour formation of torus knots in 
tightly packed DNA? Nematic arrangement in addition to 
imposing a correlation between directions of DNA 
segments that approach each other laterally also induces 
a correlation between consecutive segments of the same 
molecule, where this correlation effectively increases the 
persistence length of DNA molecules (24). To test whether 
it is the effective stiffening of DNA that causes preferential 
formation of torus knots in phage capsids, we used nu- 
merical simulations to investigate how changes of polymer 
stiffness affect knotting of modelled semiflexible polymers 
strongly confined within a sphere. 



MATERIALS AND METHODS 

Coarse-grained polymer model 

A bead-spring model based on (25) was used to model the 
polymer chain. In this model, the interaction between 
non-bonded monomers is modelled by a purely repulsive 
cut and shifted Morse potential: 

C^morse/s 

exp[-48(.T-0.8)]-2exp[-24(x-0.8)]+l if x < 0.8 
0 else 

The parameters of the model were chosen to allow com- 
parison with previous studies (26). The standard simula- 
tion unit for length is ct. At 0.8a, the interactions between 
non-adjacent beads as defined by our (repulsive) Morse 
potential become 0. Hence, we regard this value as the 
diameter of our beads. To convert to experimental units, 
we set this value equal to 2.5 nm, which is the accepted 
diameter of DNA at high salt conditions. Taking this con- 
version into account, one obtains a = 3.125 nm. 

Consecutive bonded beads are connected with FENE 
potential: 



C/pENE = -1.81n 



1 



X - 0.7 
0.3 



Using this potential, we obtain a typical distance of 
"-0.73a between consecutive bonded beads. 
Bending rigidity is introduced by: 

C/BOND(a) = BkBT{\+coi{a\), 

where a is the angle between two successive bonds. 
(In the following k^T = \ and s = 1, which corresponds 
to good solvent conditions.) A single polymer chain with 
N = 200 segments (of bead diameter 0.8a and average 
bond length of ~0.73a) is confined to a sphere of radius 
-^sphere = 5.96a. This particular value was chosen as it 
permits to have practically the same ratio between the 
diameter of the enclosing sphere and the diameter of 
modelled polymer chains as the ratio between the inner 
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diameter of phage capsids and the diameter of enclosed 
DNA (8). 

Two types of chains are investigated: The first type is 
anchored at the surface of the enclosing sphere - similar to 
viral DNA attached to the loading channel of the capsid. 
The other type is allowed to move freely within the sphere. 
In combination with local Monte Carlo moves, we apply 
the slithering snake algorithm (27) for unanchored and 
end-bridging moves for anchored chains (28) to ensure 
efficient samphng of topology. The sphere is modelled 
by an excluded volume interaction, i.e. trial configurations 
outside the sphere are rejected automatically. For each 
value of the bending rigidity B, up to 100 independent 
runs from different starting configurations were used 
in the analysis. We monitored the average unknotting 
probabihty of independent configurations over time and 
estimated correlation times between unknotted config- 
urations in single runs to ensure that the simulation runs 
are sufficiently long and the equihbrium distribution is 
reached. Error bars in Figure 2 refer to the standard 
error of the mean of the independent runs. Note that 
our computations are quite expensive and ~20 CPU 
core years were invested in this study. 

Knot identification 

Configurations are analysed by applying the foUowing 
statistical closure (29): We randomly choose two points 
on a sphere, which is much larger than the polymer, but 
has the same centre of mass. Then, these two points are 
connected by a hne with each other and one of the points 
is connected by a line with the first monomer of the chain 
while the second point is connected with the last monomer 
of the chain. The resulting closed curve is classified topo- 
logically by computing the Alexander polynomial (30). 
The whole procedure is repeated 100 times for each con- 
figuration and the knot is defined as the majority type. 
Although the Alexander polynomial does not distinguish 
between right- and left-handed forms of chiral knots this 
distinction is not required since the generic polymer model 
applied by us does not break the symmetry and therefore 
right- and left-handed forms of chiral knots are expected 
to form with the same frequencies. 

RESULTS 

To investigate how polymer stiffness affects the spectrum 
of knots resulting from cyclization of polymers in a spher- 
ical confinement, we performed Monte Carlo simulations 
that produced equihbrated statistical ensembles of long 
linear polymer chains with different stiffness but subject 
to the same spherical confinement. Our chains were 
simulated using a specific bead spring model where con- 
secutive beads in a chain could partially interpenetrate 
(31) (for more details, see 'Materials and Methods' 
section). The chains were composed of 200 beads, 
whereas the radius of the confining sphere was set to 
7.45 bead diameters. This number roughly corresponds 
to the ratio between the inner radius of P4 capsid and 
the effective diameter of DNA at high salt conditions 
(8). Assuming that the diameter of beads corresponds to 



2.5 nm and taking into account that consecutive beads 
partially interpenetrate (bond length between adjacent 
beads is roughly 0.91 bead diameters) the contour length 
of the modelled chain corresponds to ~455nm, which is 
3.5 times shorter than DNA used to study knotting in 
phage capsids (32). This shorter length has two important 
advantages: Firstly, it decreases the concentration of the 
modelled polymers inside the sphere of confinement and 
thus does not induce the isotropic/nematic transition, 
which causes an effective increase of polymer stiffness 
(24). Secondly, smaller polymer sizes aUowed us to 
explore their available configuration space in a reasonable 
simulation time. 

Our simulations investigated two types of situations. 
In the first one, the simulated polymer chains had one 
of their ends anchored at the surface of the enclosing 
sphere, which would correspond to the situation in 
wild-type phages where one end of the linear DNA stays 
in the loading channel. In the second situation, the entire 
chain could move freely within the sphere, as it is presum- 
ably the case for DNA in tailless phage mutants. 
Equilibrium configurations of chains belonging to these 
two different situations did not show systematic statistical 
differences in their shapes and in the spectrum of formed 
knots when their stiffness was the same. To characterize 
the knot spectrum that would have resulted from the cyc- 
lization of the equilibrated configurations of modelled 
polymer chains, we utilized a stochastic, unbiased 
closure method (29). The chain stiffness was controlled 
in our simulations by changing the value of a parameter 
B (see Methods section for a more precise description), 
where the higher the B value the higher the stiffness and 
where the B value of about 20 makes our model to repro- 
duce the statistical behaviour of DNA molecules. 

Cliain stiffness stimulates knotting 

Figure 1 shows overaU shapes adopted by polymers with 
the same length but different stiffness when confined 
within a small sphere of the same diameter. We can see 
that highly flexible polymers fiU the sphere of confinement 
in a rather uniform way (Figure lA). However, as the 
stiffness of polymers increases they move toward the per- 
iphery of the sphere (Figure IC and D). This result agrees 
with the intuitive expectation and an experiment that one 
could do in the lab by progressively feeding into a spher- 
ical distillation flask a long cotton string or a relatively 
stiff wire. 

A somewhat unexpected result was obtained, though, 
when we investigated how the frequency of knotting 
changes with the stiffness of modelled chains. It is 
known that in a free space knotting of polymers increases 
with the number of statistical segments in a chain (33-35). 
Therefore, if two polymers have the same physical length 
and the same effective diameter the more flexible wiU be 
more knotted than the less flexible one. However, we 
observed the opposite trend for polymers confined to a 
restricted volume. Whereas the highly flexible polymers 
were hardly knotted the frequency of knots has risen to 
>80% as the stiffness of polymers was increased (see 
insets in Figure 2). 



5132 Nucleic Acids Research, 2012, Vol. 40, No. 11 




Figure 1. Effect of increasing stiffness on the overall shapes of polymer chains having the same length and being confined within a small sphere of 
the same diameter. (A) In case of fully flexible chains (B = 0), equilibrated configurations fill the available voluine in a rather uniform way. (B). (C) 
and (D) present representative snapshots of polymer chains with increasing stiffness and having values of the parameter B amounting to 10, 20 and 
45, respectively. Notice that as chains get stiffer they progressively adopt spool-like configurations with the centre being unoccupied. Whereas the 
configuration shown in (A) is unknotted, the remaining configurations are knotted. All snapshots show unanchored chains but anchored chains have 
very similar overall appearance. 



Spectrum of knots changes with increasing stiffness 

Seeing that increasing stiffness of confined polymers leads 
to their increased tendency to form knots, we decided to 
analyse the types of formed knots. We were especially 
interested in investigating whether the increased stiffness 
can result in producing a spectrum of knots with charac- 
teristics known for the DNA knots produced in phage 
capsids. As compared to DNA knots expected to form 
in non-confined situation, DNA knots formed in phage 
capsids showed a predominance of torus knots over 
twist knots and also showed a strongly reduced frequency 
of achiral 4i knots (8). 

Figure 2 in addition to presenting the effect of polymer 
stiffness on the overall extent of knotting (see insets) also 
shows how the relative frequency of various formed knots 
changes as a function of chain stiffness. While maximally 
flexible polymers composed of 200 segments formed so 
few knots that their statistics was not reliable, polymer 
chains with a stiffness that would correspond to 4 times 
lower persistence length than this of DNA (B = 5) formed 
a spectrum of knots which is typical for knots formed in 
unconfined polymers. The twist knots 52 were roughly 
twice more frequent than torus knots 5) whereas knots 
4i were more frequent than each of five crossing knots 
(18). However, as the stiffness of modelled polymers 



increased the torus knots 5i were formed more frequently 
than the twist knots 52 and the frequency of 4i knots 
decreased below the frequencies of individual five 
crossings knots, while 3i knot remained the dominant 
knot type. Therefore, as the stiffness of modelled 
polymers confined into a sphere is increased the 
spectrum of knots progressively shifts from this character- 
istic spectrum for unconfined polymers to a spectrum 
characteristic for DNA knots formed in phage capsids. 
As shown in Figure 2A and B, non-anchored and 
anchored modelled polymers behave essentially in the 
same way with respect to frequency of formed knots. 

It needs to be said here that for computational rea- 
sons our modelled DNA chains are more than three 
times shorter than DNA molecules circularized in phage 
capsids (32). Therefore, we should not expect our 
simulated knotting spectrum to correspond to the one 
observed experimentally (32). The best we can expect is 
to reproduce the trend, i.e. observe stimulation of forma- 
tion of torus knots and observe relative suppression of 4i 
knots. 

Formed entanglements are delocalized 

In wild-type phages, the two complementary ends do not 
meet in the capsid and thus the entanglements do not 
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Figure 2. Effect of increasing stiffness on tiie probability of knotting of 
polymer chains that are about 12 times longer than the diameter of the 
confining sphere. Notice that the initial increase of stiffness greatly 
increases the overall probability of knotting (shown in insets), 
whereas a further increase of stiffness only changes the relative ratios 
between formed knots of various types. The characteristic ratio between 
different knot types observed for modelled polymers with high stiffness 
resembles the corresponding ratios in DNA knots formed in phage 
capsids (8). (A) and (B) show the knotting probabihty profiles for 
modelled polymer chains that had two or one free ends, respectively. 
Straight lines are guides to the eyes. 



become topologically fixed. However, entanglements may 
become an obstacle if formed knots could get tightened on 
the DNA and thus interfere with DNA ejection through a 
thin channel passing through phage tail (36). Although, it 
is intuitively expected that increasing stiffness should 
decrease the chance of having tight knots, we decided to 
investigate what is the effect of increasing stiffness on the 
size of the knotted portion of the chain. 

To investigate how polymer stiffness affects the spatial 
extent of knotted portions of the simulated chains, we 
concentrated on polymer chains that formed 3i knots 
upon unbiased closure and delimited then knot cores of 
individual configurations, i.e. the portion of the chain that 
is responsible for formation of 3i knot by the entire chain. 
Searching for the borders of knots' cores we followed the 
approach of References (37-39) by treating the analysed 
configurations as frozen and removing progressively 
monomers first from one end until the remaining chain 
was still forming the original knot type and then the 
chain was further trimmed from the other end using the 
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Figure 3. Increasing stiffness decompactifies knotted portions of 
confined (unanchored) polymer chains. Probability distribution func- 
tions of radii of gyration are calculated for knotted portions of 
polymer chains forming 3i knots. Values of radii of gyrations are ex- 
pressed in nm by setting the diameter of modelled chains to 2.5 nm as 
this is the accepted diameter of DNA when it is packed within phage 
capsids (21). Inset: Number of segments occupied by the knot in trefoil 
knotted chains (of total size N = 200) as a function of stiffness. For all 
B, knots are delocalized. Straight lines are guides to the eyes. 



same criteria. Subsequently, we measured radii of gyration 
of knots' cores. Figure 3 shows probability distribution 
functions of radh of gyration of 3i cores in polymer 
chains with different stiffness. It is very well visible that 
increasing stiffness increases spatial extent of knotted 
portions, even if the actual length of knotted portions is 
hardly affected as in each case knots tend to occupy large 
sections of the whole chain (inset) and are thus 
delocalized. Therefore, despite the fact that increasing 
stiffness increases the probability of knot formation by 
the confined polymers, the formed knots are in fact 
becoming less compact than knots formed on flexible 
polymers. Translating this phenomenon to the case of 
DNA packed in phage capsids, we can conclude that 
increasing stiffness decreases formation of tightened 
knotted domains and thus facihtates DNA ejection 
during phage infection. Therefore, one could expect that 
biological systems such as phages evolved in such a way 
that permits the packed DNA to increase its effective stiff- 
ness and we will discuss how this is achieved. 



DISCUSSION 

Our simulations revealed that stiffening of polymers 
confined to a small sphere naturally leads to their 
increased probability of knotting with marked predomin- 
ance of torus types of knots. This observation can be 
surprising when one takes into account that in free space 
a more stiff polymer is less Hkely to form knots than a less 
stiff polymer with the same physical length as the knotting 
probability increases with the number of statistical 
segments in unconfined polymers (33-35). 

Can we understand why more stiff polymers are more 
hkely to form knots upon confinement? In a general sense, 
to form a knot an individual polymer chain has to thread 
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Figure 4. Threading events involving neighbouring loops in coaxially spooled polymer chain lead to entanglement that can cause polymer knotting 
when the two ends anneal. (A-E) random threading events, (E-F) ends' annealing. 



through a loop formed by another portion of the same 
polymer chain. The probability of passing through 
formed loops increases with the projected area of these 
loops. In case of very flexible polymers confined to a 
sphere such as illustrated in Figure lA, the formed loops 
are usually small and therefore the probabihty of 
threading through them by other parts of the polymer is 
low. Stiff polymers confined to a sphere with the diameter 
smaller than the persistence length of the polymer natur- 
ally form very large loops with spool character (26). The 
neighbouring loops fluctuate and if one of them contains 
the free end of the polymer this end can easily thread 
through the neighbouring loop (see Figure 4A and B). 
The portions of two loops can swap their places with 
time (Figure 4C) and after some more time the free end 
can start another threading event, which could, e.g. return 
it to its original position by completing the rotation 
around the polymer chain forming the neighbouring 
loop (Figure 4D) or by reversing the half-turn rotation. 
Even in the absence of any bias favouring right- or 
left-handed winding of neighbouring loops, one should 
expect that after many such threading events the 
polymer will acquire several entanglements that will lead 
to formation of knots upon joining of polymer ends. The 
knotting equihbrium level will depend on the capacity of 
neighbouring loops to wrap around each other and this 
eventually decreases as the loops become very stiff (see 
Figure 2). A sequence of threading events that can lead 
to creation of a 3i knot is shown in Figure 4, while more 
complex sequences of threading events involving two and 
more neighbouring loops can lead to formation of more 
complex knots. 

Can we understand why stiff polymers confined into a 
sphere form preferentially 5i torus knots over 52 twist 
knots while the contrary is the case for very flexible 



polymers? At least a part of the answer is provided by 
the fact that stiff polymers within the sphere of confine- 
ment adopt trajectories without inflection points and whfle 
such an arrangement is natural for torus knots it is rather 
atypical for twist knots. Although a 52 twist knot can be 
formed by an inflection free 'toroidal' trajectory (40), this 
trajectory requires involvement of three consecutive 
coaxial loops and needs to have a complex pattern of 
interweaving involving aU three coaxial loops and 
necessitating six crossings with a specific combination of 
signs and also of under and over passages (see Figure 5). 
A 5i torus knot, in turn, in its typical toroidal form only 
needs to involve two consecutive coaxial loops and neces- 
sitates a simpler pattern of interweaving with just five 
crossings between the two loops (see Figure 5). Taking 
the above elements into account, it seems natural to 
expect that restriction of permitted configurations to 
'toroidal', inflection-free configurations entails that such 
configurations will more frequently result in formation of 
torus type 5i knots than twist type 52 knots. 

Can we understand why 4i knots form less frequently 
than 5i or 52 knots in case of configurations without the 
inflection points (37) whereas the opposite is the case for 
configurations that can populate the entire configuration 
space (17-20,41)? Similarly to the case of 52 knot the 4i 
knot belongs to the twist family of knots and its standard 
tabular representation contains two inflection points (see 
Figure 5). It is possible though to form 4i knots starting 
from a 'toroidal' inflection-free configuration (40) but this 
requires placing in a convenient proximity three consecu- 
tive coaxial loops (see Figure 5). Therefore, the overall 
arrangement preconditioning the possibility of forming a 
4i knot is quite complex and resembles the requirements 
for the formation of a 52 knot. Although formation of a 4i 
knot requires only four 'interweavings' with a specific 
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Figure 5. Coaxially spooled inflection-free configurations are more dif- 
ficult to attain for twist knots such as 4| and 52 than for torus knots 
such as 3i and 5i. 



sequence of right- and left-handed crossings and under 
and over passages and therefore in principle achieving 
this pattern may be simpler than the specific pattern 
needed to form 52 knots with six 'interweavings'. 
However, the 52 knot is chiral and exists in two enantio- 
meric forms and therefore there are two different minimal 
'interweavings' patterns that produce 52 knots while there 
is only one minimal 'interweaving' pattern that leads to 
formation of achiral 4i knots. The above explanation is 
probably not complete as we considered only minimal 
interweaving patterns leading to formation of 52 and 4| 
knots, whereas there are many other patterns that will 
permit 'toroidal' configurations to form 52 and 4i knots, 
where some wiU involve nugatory crossings etc. However, 
we believe that restricting the available configurational 
space to 'toroidal' or rather coaxial spool configurations 
is the main factor that causes that in case of DNA 
circularized within phage capsids the 4i knots are less 
frequent than 5 1 or 52 knots. 

The points discussed above help us to understand why 
equilibrated modelled stiff polymers under spherical 



confinement produce the spectrum of knots with the char- 
acteristics known for DNA knots formed in phage 
capsids. However, why do our simulations in which we 
modelled polymer chains with DNA stiffness (B = 20) 
confined to a sphere with the size corresponding to the 
size of bacteriophage capsids produce a spectrum of 
knots that was intermediate between the one known 
from phage capsids and the one modelled using very 
flexible polymers? Our simulations seem to suggest that 
the intrinsic stiffness of DNA is not sufficient to result 
in the domination of 5i knots over 52 knots and just 
suffices to make their probability of formation nearly 
equal (Figure 3). Therefore, one is led to conclude that 
the effective stiffness of DNA in phage capsids is larger 
than the stiffness of unconfined DNA. Indeed, this is most 
hkely the case. As already mentioned, while introducing our 
model, we simulated significantly shorter DNA molecules 
than phage DNA to avoid crowding that would induce 
Onsager's isotropic/nematic transition. When the DNA is 
crowded, as it is the case in phage capsids, it spontaneously 
undergoes isotropic/nematic transition leading to the ap- 
pearance of a hquid crystalhne phase. In this phase, due 
to an ordering effect involving many laterally ahgned 
DNA segments each segment is effectively stiffened 
whereas its intrinsic mechanical rigidity is unchanged (24). 
This effective stiffening of DNA most hkely explains why in 
our simulations performed in the concentration regime 
where nematic ordering was not yet induced, we had 
to stiffen our modelled polymers well above the intrinsic 
stiffness of DNA to be able to reproduce the main charac- 
teristics of the knot spectrum formed by DNA molecules 
whose ends could anneal within phage capsids. 

Our observation that effective stiffening of DNA due 
to nematic ordering is instrumental in producing the char- 
acteristic spectrum of DNA knots in phage particles is in 
agreement with a recent simulations study of Marenduzzo 
et al. (21) who also invoked a need for nematic arrange- 
ment to reproduce knotting characteristics known from 
biochemical characterization of DNA knots produced in 
phage particles. However, whereas Marenduzzo et al. (21) 
were suggesting that it is the nematic twist angle known 
for cholesteric phase of DNA liquid crystals that favours 
formation of torus knots and disfavours formation of 
achiral knots, we show here that the mere stiffening of 
DNA due to nematic arrangement is sufficient to explain 
why torus knots are favoured during DNA knot forma- 
tion in phage particles. 

In an earlier study aimed to understand the underlying 
physical principles responsible for the characteristic 
spectrum of DNA knots formed in phage particles, 
Arsuaga and Diao (42) analysed knots formed by a 
subset of equilateral random walks forming spool-like 
inflection-free configurations. The studied equilateral 
random walks had no thickness and segments were 
phantom with respect to each other. The analysis of 
formed knots in this somewhat abstract model of DNA 
packing in a phage capsids revealed though that 5i torus 
knots were just shghtly more frequent than 52 twist knot 
and that 4i knot was less frequent than 5i or 52 knots and 
therefore this simple system was close to reproducing 
the spectrum of knots formed in bacteriophage capsids. 
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The authors concluded that their model partially describes 
the situations occurring within phage capsids. We agree 
with this conclusion and show here that stiffening of DNA 
such as induced due to formation of nematic arrangement 
ensures that DNA in phage capsids adopts preferentially 
spool-like inflection-free configurations. Such a DNA ar- 
rangement, in turn, favours formation of torus type of 
knots when two ends of packed phage DNA anneal with 
each other. 

Several earlier studies proposed that low occurrence of 
achiral knots 4; among the DNA knots formed in phage 
capsids can be explained by symmetry breaking mechan- 
isms such as DNA writhe induced during DNA loading 
(8) or a preferential cholesteric twist angle between inter- 
acting DNA segments (21). We show here, however, that 
such symmetry breaking mechanisms are not needed and 
that observed suppression of achiral 4i knots can simply 
result from effective stiffening of DNA due to nematic 
ordering of tightly packed DNA. The effective stiffening, 
in turn, causes the DNA to adopt geometry of coaxial 
spools and this geometry is not favourable for formation 
of twist type of knots such as knot 4;. 

ACKNOWLEDGEMENTS 

We gratefully acknowledge computing time at the JSC 
Jiilich and the ZDV Mainz. P.V. would hke to acknow- 
ledge fruitful discussions with C. Micheletti, K. Binder 
and K. Landfester. 



FUNDING 

Swiss National Science Foundation (31003A_138367 to 
A.S.); Deutsche Forschungsgemeinschaft (SFB 625-A17 
to P.V. and D.R.) and Slovak Republic Scientific Grant 
Agency (VEGA) (2/0144/09 and 2/0093/12 to P.C.). The 
open access publication charge for this paper has been 
waived by Oxford University Press - NAR Editorial 
Board members are entitled to one free paper per year 
in recognition of their work on behalf of the journal. 

Conflict of interest statement. None declared. 
REFERENCES 

1. Wikoff.W.R., Liljas.L., Duda.R.L., Tsuruta.H., Hendrix.R.W. 
and JohnsonJ.E. (2000) Topologically linked protein rings in 
the bacteriophage HK97 capsid. Science, 289, 2129-2133. 

2. Martin.C.S., Burnett,R.M., de Haas.F.. Heinkel,R., Rutten.T., 
Fuller.S.D., Butcher,S.J. and Bamford.D.H. (2001) Combined 
EM/X-ray imaging yields a quasi-atomic model of the 
adenovirus-related bacteriophage PRDl and shows key capsid 
and inembrane interactions. Structure, 9, 917-930. 

3. Morais.M.C, Choi.K.H., Koti,J.S., Chipman,P.R., Anderson.D.L. 
and Rossniann,M.G. (2005) Conservation of the capsid structure 
in tailed dsDNA bacteriophages: the pseudoatomic structure of 
phi29. Mol Cell, 18, 149-159. 

4. Jiang,W., Chang,J., Jakana.J., Weigele.P., King,J. and Chiu,W. 
(2006) Structure of epsilonl5 bacteriophage reveals genome 
organization and DNA packaging/injection apparatus. Nature, 
439, 612-616. 

5. Comolli.L.R., Spakowitz.A.J., Siegerist,C.E., Jardine.P.J., 
Grimes.S., Anderson,D.L., Bustamante,C. and Downing,K.H. 
(2008) Three-dimensional architecture of the bacteriophage phi29 



packaged genome and elucidation of its packaging process. 
Virology, 371, 267-277. 

6. Arsuaga.J., Tan.R.K., Vazquez,M., Sumners,D.W. and 
Harvey, S.C. (2002) Investigation of viral DNA packaging 
using molecular inechanics models. Bioplivs. Cliem., 101-102, 
475^84. 

7. Petrov,A.S., Boz,M.B. and Harvey,S.C. (2007) The conformation 
of double-stranded DNA inside bacteriophages depends on capsid 
size and shape. J. Struct. Biol., 160, 241-248. 

8. Arsuaga.J., Vazquez.M., McGuirk,P., Trigueros,S., Sumners,D. 
and Roca.J. (2005) DNA knots reveal a chiral organization of 
DNA in phage capsids. Proc. Natl Acad. Sci. USA, 102, 
9165-9169. 

9. Yu,J., Moffitt,J., Hetherington,C.L., Bustamante,C. and Oster,G. 
(2010) Mechanochemistry of a viral DNA packaging motor. 

/. Mol. Biol.. 400, 186-203. 

10. Kornyshev,A.A. and Leikin.S. (2000) Twist in chiral interaction 
between biological helices. Fhy.i. Rev. Lett., 84, 2537-2540. 

11. Chattoraj,D.K. and Inman,R.B. (1974) Location of DNA ends in 
P2, 186, P4 and lambda bacteriophage heads. /. Mol. Biol., 87, 

1 1-22. 

12. Wang,J.C., Martin,K.V. and Calendar,R. (1973) On the sequence 
similarity of the cohesive ends of coliphage P4, P2, and 186 
deoxyribonucleic acid. Biocheini.itry, 12, 2119-2123. 

13. Liu.L.F., Perkocha.L., Calendar,R. and Wang,J.C. (1981) 
Knotted DNA from bacteriophage capsids. Proc. Natl Acad. Sci. 
USA, 78, 5498-5502. 

14. Vologodskii,A.V., Crisona,N.J., Laurie,B., Pieranski,P., 
Katritch.V., Dubochet,J. and Stasiak.A. (1998) Sedimentation 
and electrophoretic migration of DNA knots and catenanes. 
J. Mol. Biol, 278, 1-3. 

15. Trigueros.S., Arsuaga.J., Vazquez.M. E., Sumners.D.W. and 
Roca.J. (2001) Novel display of knotted DNA molecules by 
two-dimensional gel electrophoresis. Nucleic Acids Res., 29, 
E67-E67. 

16. Arsuaga,J., Vazquez.M., Trigueros.S., Sumners.D. and Roca.J. 
(2002) Knotting probability of DNA molecules confined in 
restricted volumes: DNA knotting in phage capsids. Proc. Natl 
Acad Sci. USA. 99. 5373-5377. 

17. Tsurusaki.K. and Deguchi.T. (1995) Fractions of particular 
knots in Gaussian random polygons. J. Phv.s. Soc. Jpn., 64. 
1506-1518. 

18. Katritch.V.. Olson.W.K.. Vologodskii.A.. Dubochet.J. and 
Stasiak.A. (2000) Tightness of random knotting. Phy.s. Rev. E. 
Stat. PhY.s., 61. 5545-5549. 

19. Shaw.S.Y. and Wang,J.C. (1993) Knotting of a DNA chain 
during ring closure. Science, 260, 533-536. 

20. Rybenkov,V.V., Cozzareni,N.R. and Vologodskii.A.V. (1993) 
Probabihty of DNA knotting and the effective diameter of the 
DNA double helix. Proc. Natl Acad. Sci. USA, 90, 5307-5311. 

21. Marenduzzo.D., Orlandim,E., Stasiak.A.. Sumners de.W., 
Tubiana.L. and Miclieletti,C. (2009) DNA-DNA interactions in 
bacteriophage capsids are responsible for the observed DNA 
knotting. Proc. Natl Acad. Sci. USA, 106, 22269-22274. 

22. Marenduzzo.D., Micheletti, C. and Orlandini.E. (2010) Biopolymer 
organization upon confinement. /. Phv.s. Condens. Matter. 11, 
283102. 

23. Stanley,C.B., Hong,H. and Strey,H.H. (2005) DNA cholesteric 
pitch as a function of density and ionic strength. Biophr.i. J.. 89, 
2552-2557. 

24. Lin.S.. Numasawa.N.. Nose.T. and Lin.J. (2007) Corse-grained 
molecular dynamic simulations for lyotropic liquid-crystalline 
solutions of semiflexible rod-like molecules. Mol. Cry.it. Liq. 
Cry.st., 466, 53-76. 

25. Cifra,P., Benkova,Z. and Bleha.T. (2008) Persistence lengths and 
structure factors of wormlike polymers under confinement. 

/. Phys. Cheni. B, 112, 1367-1375. 

26. Cifra,P. and Bleha.T. Shape transition of semi-flexible 
macroniolecules confined in channel and cavity. Eur. Phys. J. E., 
32. 273-279. 

27. Kron.A.K.. Ptitsyn.O.B.. Skvortsov.A.M. and Federov.A.K. 
(1967) A study of statistical globula-coil transition in 
macromolecules using the Monte-Carlo technique. Molek. Biol., 1. 
576-582. 



Nucleic Acids Research, 2012, Vol. 40, No. 11 5137 



28. Reith,D. and Virnau,P. (2010) Implementation and performance 
analysis of bridging Monte Carlo moves for off-lattice single 
chain polymers in globular states. Coinput. Phys. Commun., 181, 
800-805. 

29. Mansfield,M.L. (1994) Are there knots in proteins? Nat. Struct. 
Biol., 1, 213-214. 

30. Virnau,P. (2010) Detection and visualization of physical knots in 
macromolecules. Phy.i. Procedia, 6, 117-125. 

31. Cifra,P., Benkova.Z. and Bleha,T. (2008) Persistence lengths and 
structure factors of wormlike polymers under confinement. 

/. Plm. Chem. B, 112, 1367-1375. 

32. Trigueros,S. and Roca,J. (2007) Production of highly knotted 
DNA by means of cosmid circularization inside phage capsids. 
BMC Biotechnol., 7, 94. 

33. Frisch,H.L. and Wasserman,E. (1961) Chemical topology. 
/. Am. Chem. Soc, 83, 3789-3795. 

34. Delbruck,M. (1962) Mathematical Problem.^ In the Biological 
Science.^. Vol. 14. Mathematical Society, Providence, Rhode 
Island, pp. 55-63. 

35. Frank-Kamenetskii,M.D., Lukashin,A.V. and Vologodskii,A.V. 
(1975) Statistical mechanics and topology of polymer chains. 
Nature, 258, 398^02. 



36. Matthews,R., Louis,A.A. and Yeomans,J.M. (2009) 
Knot-controlled ejection of a polymer from a virus capsid. Phy.s. 
Rev. Lett., 102, 088101. 

37. Marcone,B., Orlandini,E., Stella,A.L. and Zonta,F. (2007) 
Size of knots in ring polymers. Phjs. Rer. E Stat. Nonlin. Soft 
Matter Phys., 15, 041105. 

38. Virnau.P., Kantor,Y. and Kardar,M. (2005) Knots in globule and 
coil phases of a model polyethylene. /. Am. Chem. Soc, 127, 
15102-15106. 

39. Virnau,P., Mirny, L.A. and Kardar,M. (2006) Intricate knots 
in proteins: Function and evolution. PLoS Comput. Biol., 2. 
el22. 

40. Ricca,R.L. (1998) New developments in topological fluid mechanics: 
from Kelvin's vortex knots to magnetic knots. In: Stasiak,A., 
Katritch,V. and Kauffman,L.H. (eds). Ideal Knots, Vol. 19. World 
Scientific, Singapore, pp. 255-273. 

41. Micheletti,C., Marenduzzo,D., Orlandini,E. and Sumners,D.W. 
(2008) Simulations of knotting in confined circular DNA. 
Biophys. J., 95, 3591-3599. 

42. Arsuaga,J. and Diao,Y. (2008) DNA knotting in spooling like 
configurations in bacteriophages. Comp. Math. Methods Med., 9, 
303-316. 



